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1. Introduction 

The search for the topology of the spatial sections of the universe has made tremendous 
progress in the past years (see Lachieze-Rey and Luminet (1995) for an introduction to 
the subject and early references, Luminet and Roukema (1999) and Uzan et al. (1999b) 
for a review of the late developments). Methods using two-dimensional data sets, such as 
the cosmic microwave background maps planned to be obtained by the MAP and Planck 
Send offprint requests to: Roland Lehoucq 
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Surveyor satellite missions, and three-dimensional (3D) data sets, such as galaxy, cluster 
and quasar surveys with redshifts, have been developped. 

Following our previous works (Lehoucq et al. 1996, Lehoucq et al. 1999 and Uzan 
et al. 1999a), we focus on the topological information that can be extracted from a 
3D catalogue of cosmic objects. The key point of all the topology-detecting methods 
is based on the "topological lens effect", i.e. on the fact that if the spatial section of 
the universe has at least one characteristic size smaller than the spatial scale of the 
catalogue, then one should observe different images of the same object. The original 
crystallographic method (Lehoucq et al. 1996) used a pair separation histogram (PSH) 
depicting the number of pairs of catalogue's objects having the same three-dimensional 
separations in the universal covering space, with the idea that spikes should stand out 
dramatically at characteristic lengths related to the size of the fundamental domain and 
to the holonomies of space. However, following a remark by Weeks (1998), we proved 
(Uzan et al. 1999a and Lehoucq et al. 1999) that sharp spikes can emerge in the PSH 
only if the holonomies of space are Clifford translations - a result independently derived 
by Gomero et al. (1998). As a consequence, the PSH method docs not apply to the 
detection of topology in universes with hyperbolic spatial sections. 

Since then, various generalisations of the crystallographic method were proposed. Fa- 
gundes and Gaussman (1998) suggested to map the differences between the PSH of a 
simulated catalogue in a compact hyperbolic universe and the PSH in the correspond- 
ing simply-connected universe having the same distribution of objects and cosmological 
parameters. They noticed sharp oscillations on the scale of the bin width, modulated by 
a broad oscillation on the scale of the curvature radius. However in Uzan et al. (1999b), 
we calculated on one hand the differential PSH between a simulated distribution in 
a compact hyperbolic Weeks space and a similar distribution in the simply-connected 
hyperbolic space H 3 , on the other hand the differential PSH between two different dis- 
tributions (with the same total number of objects) in the simply-connected hyperbolic 
space H 3 . Both curves (fig. 10 of Uzan et al. 1999b) exhibit the same pattern of sharp 
and broad oscillations, which shows that the topological significance of such a pattern is 
highly doubtful. 

Fagundes and Gaussman (1999) next proposed a modified crystallographic method in 
which the topological images in simulated catalogues are pulled back to the fundamental 
domain before the set of 3D distances is calculated. The distribution of pair distances 
is expected to be peaked around zero. The main drawback of this method is that the 
nature of the signal strongly depends on the topological type, on the orientation of the 
fundamental domain and on the position of the observer inside the latter. Thus the 
pullback method can be useful only if the exact topology is already known. On the other 
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hand, Gomero et al. (1999) introduced a mean PSH, aimed to reduce the statistical noise 
that may mask the topological signature. They claimed that such a technique should 
detect the contribution of non-translational isometries to the topological signal, whatever 
the curvature of space, but the applicability of their method to real data has not been 
demonstrated. 

In order to improve the signal-to-noise ratio inherent to PSH's , instead of reduc- 
ing the statistical noise we can enhance the topological signal by collecting all distance 
correlations into a single index. Such was our purpose in Uzan et al. 1999a, where we 
reformulated cosmic crystallography as a collecting-correlatcd-pairs method (CCP). The 
CCP technique rests on the basic fact that in a multiply connected universe, equal dis- 
tances appear more often than just by chance, whatever the curvature of space and the 
nature of holonomies. We also showed that the extraction of a topological signal drasti- 
cally depends on the rather accurate knowledge of the cosmological parameters, namely 
the density parameter, O , and the cosmological constant parameter, Qao- This is due to 
the fact that all cosmic crystallography methods require the determination of 3D separa- 
tions between two any images: observations use redshifts for determining the coordinate 
distance in the universal covering space (in addition to the angular positions on the ce- 
lestial sphere), and the redshift-distance relation involves the cosmological parameters. 
Conversely, the detection of a topological signal would help to determine accurately the 
curvature parameters (see also Roukema and Luminet 1999). 

However it is to be recognized that, in the framework of cosmic crystallography, 
the application of numerical simulations to real data rest on two idealized assumptions, 
namely: 

1. all objects are strictly comoving 

2. the catalogues of observed objects (quasars, galaxy clusters) are complete. 

Although quoted in Lachieze-Rey and Luminet (1995), the quantitative effects of 
such simplifications have not been fully discussed in the literature. In Lehoucq et al. 
(1996), the authors qualitatively argued that the distortion due to peculiar velocities 
was negligible, and they performed numerical simulations to study the influence of the 
angular resolution of the surveys. In Roukema (1996), the influence of the astrophysical 
uncertainties (mainly of spectroscopic measurements and of peculiar velocities) on a 
method trying to find quintuplets of quasars with the same geometry (and thus which 
may be topological images) was evaluated. It was shown that, in such a case, the most 
serious uncertainty comes from the radial peculiar velocities. In Uzan et al. (1999a), we 
discussed the effect of the peculiar velocities and of the errors in the spatial position 
arising from the imprecisions on the values of the cosmological parameters. To finish, the 
errors on the determination of the position and of the peculiar velocities were discussed 
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in Roukema (1996), and the way that the constraints on f2o and Oao depend on the 
redshifts of multiple topological images and on their radial and tangential separations 
was calculated in Roukema and Luminet (1999). 

The goal of the present article is to have a critical attitude on the methods we have 
developped so far, by listing all the sources of observational uncertainties and by eval- 
uating their effects on the theoretical efficiencies of the PSH method (Lehoucq et al. 
1996) and of the CCP method (Uzan et al. 1999a) respectively. In § ||, we discuss the 
nature of uncertainties and in § ^ the numerical methods used to estimate their effects 
on the topological signal. We then quantify the magnitude of each effect in the Euclidean 
and hyperbolic cases (§ ||) (we postpone the case of elliptic spaces to a further study). 
In conclusion (§ ||) we compare our results to the performances of current and future 
observational programs aimed to detect the topology of the universe. 



Notations and descriptions 

We keep the notations of our previous articles (Lehoucq et al. 1999, Uzan et al. 1999a). 
The local geometry of the universe is described by a Friedmann-Lemaitre metric 

ds 2 = -dt 2 + a 2 (t) [d X 2 + f(x) (di? 2 + sin 2 tfd^ 2 )] , (1) 
where a is the scale factor, t the cosmic time, x the comoving radial distance, and 
fix) — ( sm X> X; smn x) according to the sign of the space curvature k = (+1, 0, —1). The 
time evolution of a is obtained by solving the Friedmann and the conservation equations 

3 a z 3 

p = -3H(p + P) (3) 

where k = 8irG/c 4 , p and P are respectively the matter density and the pressure, A is 
the cosmological constant and H = a/a is the Hubble parameter (with a dot referring to 
a time derivative). We also use the standard parameters 

™ d (4 » 

This completely specifies the properties and the dynamics of the universal covering space. 
In the following, we assume that we are in the matter dominated era, so that P = and 
p oc a~ 3 . 

The topology of the spatial sections is described by the fundamental domain, a poly- 
hedron whose faces are pairwise identified by the elements g of the holonomy group T 
(see Lachieze-Rey and Luminet (1995) for the details). 



2. Observational uncertainties 

Early methods used to put a lower bound on the size of the universe and based on the 
direct recognition of multiple images of given objects - such as our Galaxy (Sokolov 
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and Schvartsman 1974), the Coma cluster (Gott 1980) - had to face a major drawback 
due to the fact that the same object would be seen at different lookback times. Thus, 
evolution effects such as photometric or morphologic changes will, in most cases, render 
the identification of objects impossible (see however Roukema and Edge, 1997). Statistical 
techniques such as the PSH method (Lehoucq et al. 1999) and the CCP method (Uzan 
et al. 1999a) are free from such biases. 

The two main sources of uncertainties for all the statistical methods trying to detect 
the topology of the universe in catalogues of cosmic objects are 

(A) the errors in the positions of observed objects, which can be separated into: 

(Al) the uncertainty in the determination of the redshifts due to spectroscopic impreci- 
sion; such an effect is purely experimental and exists even if the objects are strictly 



(A2) the uncertainty in the position due to peculiar velocities of objects, which induce 

peculiar redshift corrections 
(A3) the uncertainty in the cosmological parameters, which induces an error in the 

determination of the radial distance (via the redshift - distance relation) 
(A4) the angular displacement due to gravitational lensing by large scale structure. 
(B) the incompleteness of the catalogue, which has two main origins: 

(Bl) selection effects implying that some objects are missing from the catalogue 
(B2) the partial coverage of the celestial sphere, due either to the presence of the galactic 

plane, or to the fact that surveys are performed within solid angles much less than 



Effects of peculiar velocities 

A peculiar velocity has two effects on the determination of the position of a cosmic 
object: 

(i) An integrated effect, coming from the fact that if a galaxy has a proper velocity, then 
its true position differs from its comoving one. 

(ii) An instantaneous effect, due to the fact that the radial component of the proper 
velocity will add to the cosmological redshift an extra term, Azd op - 

Assuming that the vector velocity v of a galaxy is constant, which is indeed a good 
enough approximation to estimate the effects of peculiar velocities, a galaxy at redshift 
z has moved from its comoving position by a comoving distance 



where t[z] is the look-back time, obtained by integrating the photon geodesic equation 



comoving 



47T. 



51 = vt[z] 



(5) 




(6) 



6 
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Fig. 1. Variation of Az pos in units of («||/c) as a function of the redshift, assuming the 
galaxy peculiar velocity is constant. 



This reduces to the well known expression 
2 (. 1 



1 - 



(7) 



3H y (1 + z)3/ 2 
when fio = 1 and ^ao = 0. 

Now, 6£ can be expressed in terms of a peculiar redshift, A2p OS , and of an angular 

displacement, A9 pos , which depend on the galaxy velocity as 

f || ct[z] 



A 



AO- 



c x [z\ 
\v± \ ct[z] 



(8) 

where a prime refers to derivative with respect to z, Vu is the component of the velocity v 
with respect to the line-of-sight direction 7, and v± is the non-radial peculiar velocity, 
defined as 



v.-f and Vj 



(10) 



(11) 



x[z] is the observer area distance, given by 
x[z] = -2-? ^ 

a H J_i_ y/n A0 x 4 + (1 - n - n A0 )x 2 + n x 

In figures [l] and m, we respectively depict the variations of Az pos and A6> pos as a function 
of the redshift. 

Concerning the instantaneous effect ii) , the redshift uncertainty AzDop can be related 
to the galaxy proper velocity as follows. 

If we consider the trajectory of a photon x fl (s), s being the afiine parameter along 
the null geodesic, the relation between the emission (E) wavelength Ae and the reception 
(R) wavelength Ar can be expressed as 

Ae 



(12) 
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Fig. 2. Variation of A# pos as a function of the redshift. As long as z < 5, we have that 
Af? pos /lrad < 0.5(v±/c). Note also that when z — > 0, A0 pos /lrad — ► (vx/c). 



where fc M = dx^/ds is the tangent vector to the photon geodesic; it° bs and ujf 1 are 
respectively the 4-velocity of the observer and of the galaxy. Neglecting the perturbations 
of the metric and of the matter and focusing on the Doppler effect, one can easily show 
that 

fc"« M = — ^[i + V-l, (13) 
a(t) L c J 

where 7* is the direction in which the galaxy is observed and v l is its (Newtonian) 
velocity. The observed (spectroscopic) redshift, z b s , and the cosmological redshift, z cosm , 
respectively defined by 

_ a(t E ) 



11 — ^ R 11 
1 + z ohs = — , 1 + / 

are thus related by 

1 ~\~ Zq]j S (1 -\- 2 cosm ) 



1 + r 



a(tn) 



,,obs\ 



(14) 



(15) 



Assuming that we can substract the Earth's velocity, an error Az = z Q b s — z cosm in the 
determination of z D b s can be interpreted in terms of a radial peculiar velocity given by 



Azoop = ^-(l + z c 
c 



(16) 



As expected, a galaxy receding from us, i.e. with 7 I i>f al > will have an additional red- 
shift, whereas a galaxy drawing nearer to us will induce a Doppler blueshift correction. 



The considerations above are also useful for discussing the effects of angular dis- 
placement A6* g i due to gravitational lensing by large scale structure (Mellier 1999, van 
Waerbeke et al. 2000). A typical value A8 g i ~ 1 arcsec corresponds to a peculiar velocity 
a few km.s -1 as long as the redshift is less than 5 (as can be seen from Fig. |^). Thus 
effect (^44) is much less important than effect (A2). 
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Effects of catalogue incompleteness 

Objects can be missing from a catalogue due to strong evolution effects (this is par- 
ticularly the case for quasars), and to absorption of light by intergalactic gas or dust in 
some sky directions. Another well known selection effect is the Malmquist bias: statisti- 
cal samples of astronomical objects which are limited by apparent magnitude have mean 
absolute magnitudes which are different than those of distance-limited samples. An ap- 
parent magnitude-limited sample contains, if luminosity function has a finite width, some 
very luminous objects which, in spite of their large distances, can jump the apparent- 
magnitude limit of the catalogue. While one looks further and further out one finds more 
and more luminous objects. The consequence for an apparent-magnitude catalog (which 
is the case with galaxy catalogues) is that dwarf galaxies fade out quickly with distance, 
and finally at the largest distances the extremely luminous and rare galaxies are the only 
ones which can enter the catalogue. 

In addition to effects (Bl) and (B2), some noise spikes may appear due to gravita- 
tional clustering of objects. The large scale distribution of galaxies shows a variety of 
landscapes containing voids, walls, filaments and clusters in a complex 3D sponge-like 
pattern. For instance, since there are many galaxies in clusters, the distances associated 
to cluster-cluster separations may appear as fake spikes in the PSH for galaxies. Typical 
cluster - cluster separations are lSOh^ 1 Mpc (Guzzo et al. 1999). Such effects currently 
occur in N-body simulations that show clustering (Park and Gott 1991). The trouble 
can be relieved if galaxy clusters instead of galaxies are used as typical objects for prob- 
ing the topology, although in that case supercluster-supercluster separations might also 
introduce noise spikes (at a lower level). 



3. Numerical implementation 

In the following, we perform numerical simulations to evaluate separately the magnitudes 
of the effects listed above. As usual, we start from a random distribution of cosmic objects 
in the fundamental domain. In a first step we generate a complete catalogue of comoving 
objects by unfolding the distribution in the universal covering space. We refer to this 
catalogue as the ideal catalogue since it would correspond to ideal observations. In a 
second step we introduce the various errors {A{) — (Bj) in order to build more realistic 
catalogues which depart from the ideal one. 
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3.1. Uncertainties on positions 

To study numerically the errors on the positions, we first assume that the cosmological 
parameters are known with good enough accuracy, so that we do not discuss (A3). Then 
we perform the following calculations. 

1. To evaluate the effect of the observational imprecision, we give to each object of the 
ideal catalogue (built on the assumption that the objects are strictly comoving) a 
redshift error to be added to its ideal redshift. We assume that the distribution of the 
redshift error is Gaussian, with mean value z = and dispersion Az. 

Note that, since Az is absolute, the relative error will be less important when we deal 
with catalogues of higher redshift objects. 

2. To discuss (A2), we associate a peculiar velocity to each point of the catalogue before 
unfolding. Hence, there will be correlations between the velocities of two topological 
images. It follows that, assuming that the time evolution of the peculiar velocity is 
small, the observed velocities of two topological images of the same object can have 
different directions but similar magnitudes (see e.g. Roukema and Bajtlik 1999). As 
in the previous case, we assume that the velocity distribution is Gaussian. Moreover, 
velocities at different points of space are assumed to be uncorrelated both in magni- 
tude and direction, which is indeed not strictly the case in real data since large scale 
streaming motions of galaxies have been observed (Strauss and Willick 1995). 

More precisely, to generate the catalog of topological images taking into account the 
peculiar velocities, we proceed as follows: 

1. We generate a random collection of points Mi (named original sources) uniformly dis- 
tributed inside the fundamental domain. Each point is assigned a velocity v% according 
to a Gaussian distribution with mean v and dispersion Av. 

2. We unfold this catalogue of original sources to obtain the set of points Mk.i, images 
of Mi by the holonomies gt G T. 

3. Each image Mk,i has a redshift Zk,i corresponding to a look-back time = r(zk,i) 
calculated with formula @. 

4. The final catalogue accounting for the peculiar velocities of the original sources is 
obtained by applying the gk to Mi + Vi x Tk,i for all couples (k,i). 

Since Az can be interpreted either as an uncertainty on the redshift determination or 
as due to a radial peculiar velocity, both errors (Al) and (A2) can be investigated with 
the same calculations. The interpretation of Az in terms of a radial peculiar velocity is 
useful to compare the orders of magnitude of (Al) and (A2). In real data, the error due 
to (A2) is expected to be greater than the error due to (Al). 



10 R. Lehoucq, J-P. Uzan & J-P. Luminet: Limits of Crystallographic Methods 

Table 1. The multiplication factor, i.e. the ratio between the total number of entries 
in the catalogue and the number of "original" objects randomly distributed in the fun- 
damental domain, is given as a function of the solid angle of the survey. We have also 
indicated the corresponding sky coverage. 
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3.2. Catalogue incompleteness 

Starting from an ideal catalogue, we simulate two kinds of incomplete catalogues: 

1. To account for effect (Bl), we randomly throw out p% of the objects from the ideal 
catalogue. In our various runs we vary p while keeping approximately constant the 
number of catalogue objects, which means that we must increase the number of 
"original" objects in the fundamental domain when p is increasing. 

2. To simulate effect (B2), we generate a catalogue limited in solid angle by selecting 
only the objects that lie in a beam of aperture 9. The sky coverage is related to 9 via 

q% = \ {l- cos(0/2)) . (17) 
Again, we keep constant the number of objects in the catalogue when we vary 9. 
Table [I] below gives typical numbers. 

As a matter of fact, more than the aperture angle, the depth of the survey will be 
critical for the multiplication factor, since within a beam of given angle, if the redshift 
cut-off is great enough to encompass a distance (in the universal covering space) N times 
greater than the size of the fundamental domain, at least N topological images will be 
expected in the beam. 

4. General Results 

4-1. Euclidean spaces 

We first consider the Euclidean case and we apply the PSH crystallographic method 
as described in Lehoucq et al. (1996), restricting to a typical situation where f^o = 
0.3 and f^Ao = 0.7. We choose the topology of the universe to be a cubic 3-torus T\ 
(see e.g. Lachieze-Rey and Luminet (1995) for description), with identification length 
L = 3, 000 Mpc for a Hubble parameter Hq = 75km.s -1 Mpc~ 1 . In the various runs, the 
number of objects in the catalogue is kept constant at 8500 (this is the order of magnitude 
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Fig. 3. Plot of Azi as a function of the depth of the catalogue, for a cubic hypertorus of 
size 3,000 Mpc which corresponds to 45% of the Hubble radius (4, 950 h^ 1 Mpc) when 
Oo = 0.3, f^Ao = 0.7 and h = 0.75, using the PSH method. 

Table 2. Values of the critical percentage of rejection above which the PSH spikes 
disappear, as a function of the redshift cut-off. 



%cut 
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> 3 


Pll%] 


70 
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> 90 



of the number of objects in current quasar catalogues). We examine separately the effects 
of errors in position due to redshift uncertainty Az and peculiar velocities Av, and the 
effects of catalogue incompleteness due to selection effects and partial sky coverage. Each 
of these effects will contribute to spoil the sharpness of the topological signal. For a given 
depth of the catalogue, namely a redshift cut-off z cut , we perform the runs to look for 
the critical value of the error at which the topological signal fades out. 

Figure || gives the critical redshift error Azi above which the topological spikes dis- 
appear. 

The effect of peculiar velocities is very weak, since we find that the peculiar velocity 
must exceed Avi = 10, 000 km. s~ when z cut — 1 and Avi — 40, 000 km. s - when z cut — 5 
in order to make the topological signal disappear. As already pointed out in Lachieze- 
Rey and Luminet (1995), the observed peculiar velocities of galaxies have typical values 
much less that Avi. 

The effects of catalogue incompleteness are summarized in tables || and || The topo- 
logical signal would be destroyed only for a very large rejection percentage or a small 
aperture angle. 
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Table 3. Values of the aperture angle below which the PSH spikes disappear, as a 
function of the rcdshift cut-off. 
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4-2. Hyperbolic spaces 

We now turn to universes with hyperbolic spatial sections and apply the CCP method as 
described in Uzan et al. (1999a). The obtention of a topological signal (the so-called CCP 
index) strongly depends on the correct determination of the cosmological parameters. In 
Uzan et al. (1999a) we discussed the problems arising from the spanning of the cosmolog- 
ical parameters space with a required "accuracy bin" e. Latest independent constraints 
on these cosmological parameters from the cosmic microwave background (de Bernardis 
et al. 2000), the study of supernovae (Efstathiou et al. 1999), of large scale structure at 
z = 2 (Roukema and Mamon, 2000) and of gravitational lensing (Mellier 1999), make 
us hope that we can restrict further the parameters space to apply efficiently the CCP 
method. In the following, we assume that the cosmological parameters are known with 
the required accuracy. 

We choose the topology of the universe to be described by a Weeks manifold (Weeks, 
1985) (see also Lehoucq et al. (1999) for the numerical implementation of this topology), 
assuming the cosmological parameters given by £Iq = 0.3 and f^Ao = 0. In such a case the 
number of copies of the fundamental domain within the horizon is about 190. However 
our simulated catalogues are much smaller than the horizon volume. We fix the number 
of objects to 1300, which is a compromise between a realistic catalogue population and 
a reasonable computing time. 

Again we perform the tests by varying the errors Az and Av around their average 
values z = and v = until when the CCP index falls down to noise level. 

Contrarily to the Euclidean case, both slight changes in redshift and in peculiar 
velocity induce errors in position which dramatically eradicate the topological signal: an 
error of only ~ 50kms~ in velocity and an error in redshift of the order of the bin 
accuracy e ~ 10~ 6 drown the CCP index into noise. 

The incompleteness effects are less dramatic, as shown in figures | and |. A gain, for 
a given redshift cut-off, we performed the runs by varying p and 9 in order to find the 
critical values at which the topological signal disappears. 
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Fig. 4. Plot of the CCP-index (normalized to background noise) as a function of the 
rejection percentage of objects pi, for various values of of the depth of the catalogue. The 
topological signal disappears when the rejection percentage is greater than pi . 
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Fig. 5. Plot of the CCP-index as a function of 9i for various values of the depth of the 
catalogue. The topological signal disappears when the aperture angle falls down below 

0i. 



5. Conclusions and perspectives 



Our numerical results have now to be compared with the precisions of present experi- 
mental 3D data, and to the performances of observational programs started or expected 
to be achieved in the next decade. 

At present day, a typical precision practical for the spectroscopic uncertainty is Az ~ 
0.001 for an object such as a quasar. In clusters, spectroscopic redshifts can be found 
very precisely for individual galaxies. 
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Concerning peculiar velocities, the typical dispersion velocity is 1000 km. s -1 in rich 
clusters. The X-ray velocity of the peak of the X-ray distribution would also provide a 
way to estimate the true cluster redshift, including the peculiar velocity of the cluster 
as a whole. For a quasar, a conservative upper limit to the peculiar velocity, assuming 
the quasar to be at the centre of a galaxy, can be taken as Az ~ 0.002. So, from an 
experimental point of view, the uncertainties on the redshifts will be dominated by effect 
(A2), namely the peculiar velocities, rather than by the spectroscopic imprecision. 

The main limitation of present 3D samples is the small volume of existing redshift 
data. Future surveys will significantly improve both the redshift cut-off and the sky 
coverage. For instance, the Sloan Digital Sky Survey (SDSS) (Loveday 1998) will map in 
detail one-quarter of the entire sky, determining the positions and absolute brightnesses 
of more than 100 million celestial objects. It will also measure distances to more than a 
million galaxies and quasars. More precisely, SDSS will map a contiguous ir steradians 
area in the north Galactic cap, up to a limiting magnitude 23 for two thirds of the 
observing time, together with three southern stripes centred at RA a — 5° and with 
central declinations of 5 = +15°, 0° and —10° for the remaining one third of the time. 
Concerning the distance determinations, 10 6 galaxies and 10 5 quasars will be observed 
spectroscopically with a resolution Az/z ~ 5.10~ 4 . The main galaxy sample will consist 
of ~ 900,000 galaxies up to magnitude 18, with a median redshift z ~ 0.1. A second 
galaxy sample will consist in ~ 100, 000 luminous red galaxies to magnitude 19.5 with a 
median redshift z ~ 0.5. Precision on redshifts can be estimated for the reddest galaxies 
to Az <~ 0.02. A sample of ~ 100, 000 quasars will be observed, an order of magnitude 
larger than any existing quasar catalogue. The complete survey data will become public 
by 2005. We can also mention the ESO-VLT Virmos Deep Survey Project (Lefevre 2000), 
a comprehensive imaging and redshift survey of the deep universe based on more than 
150 000 redshifts. 

In the field of X-ray observations, the XMM satellite (Arnaud 1996) will provide deep 
insight on X-ray galaxy clusters and active galactic nuclei. Also, the XEUS project under 
study by ESA (Parmar et al. 1999) will be a long-term X-ray observatory at 1 keV, with 
a limiting sensitivity around 250 times better than XMM, allowing XEUS to study the 
properties of galaxy groups at z = 2 and active galactic nuclei at z ~ 3. 

In the present paper we have investigated how the various observational uncertainties 
will spoil the topological signal expected to arise in the ideal situation when crystallo- 
graphic methods are applied to complete catalogues of perfectly comoving objets with 
zero proper velocities and whose 3D-positions are known with infinite accuracy. By nu- 
merical simulations we have introduced random errors for each possible uncertainty, and 
we varied the parameters to determine the limits at which the topological signal vanishes. 
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Our numerical calculations of the spoiling effects due to the various uncertainties 
(Ai) — (Bi) clearly show that the crystallographic methods are stable (in the sense that 
the topological signal is robust when data depart from the ideal ones) in the Euclidean 
case, but highly unstable in the hyperbolic case. Indeed in a small multi-connected flat 
space, realistic values of peculiar velocities of objects, errors in rcdshift determinations 
and partial sky coverage will not make the PSH method to fail. This can be understood 
by the fact that the topological images of a given object are related together by Clifford 
translations which enhance the topological signal. On the contrary, in a compact hyper- 
bolic model, holonomies are not Clifford translations. The topological signal, built as a 
CCP index, is destroyed as soon as small errors are introduced in the position of objects, 
due either to peculiar velocities or to redshift measurement imprecision. 

The same kind of critical analysis should be made with the 2D topology-detecting 
methods based on the analysis of CMB data. For instance, in the pairs of matched circles 
method (Cornish et al. 1998), it would be necessary to investigate how deviations to the 
ideal situation, such as a non zero thickness of the last scattering surface or peculiar 
motions of the emitting primordial plasma regions, would alter the pattern of perfectly 
matched circles. 

To conclude, let us comment on the future of observational cosmic topology, at 
the light of observational constraints on the curvature parameters recently provided by 
BOOMERanG and MAXIMA balloon measurements (de Bernardis et al. 2000, Hanany et 
al. 2000). Under specific assumptions such as a cold dark matter model and a primordial 
density fluctuation power spectrum, the range of values allowed for the energy-density 
parameter £1 = Sl + ^ao is restricted to 0.88 < £1 < 1.12 with 95% confidence. This 
means that the curvature radius of space is as least as great as the radius of the observ- 
able universe (delimited by the last scattering surface). Such results, if accepted, leave 
open all three cases of space curvature as well as most of multi-connected topologies. 

A strictly flat space is quite improbable. Even inflationary scenarios predict a value 
of £1 asymptocally close to 1, but not strictly equal. From a topological point of view, a 
strictly Euclidean space would be interesting since we have shown that the PSH method 
is robust enough to provide a topological signal even when realistic uncertainties on the 
data are taken into account. 

For compact hyperbolic spaces, if we accept the recent observational constraints 
0.88 < Oo + Oao < 1, a number of spaceforms such as the Weeks or the Thurston 
manifolds still have topological lengths smaller than the horizon size (Weeks: SnapPea). 
However the topological lens effects would be weaker in spaceforms with £1 close to 1 
than in spaceforms with £1 ~ 0.3 (see figs. 2 and 3 of Lehoucq et al. 1999). Furthermore, 
only the CCP method can be applied for detecting the topology, and the present article 
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shows that the topological signal will fall to noise level as soon as uncertainties smaller 
than the experimental ones are taken into account in the simulations. 

Eventually, elliptical spaceforms appear to be the most interesting case. On a theo- 
retical point of view, as far as we know, no inflationary model is able to drive the density 
parameter to a value greater than 1, so that if space happened to be really elliptical, new 
models should be built in order to explain, e.g., the primordial fluctuations spectrum or 
the horizon problem. From a topological point of view, since the volumes of (all closed) 
spaceforms are not bounded below (the order of the holonomy group can be arbitrarily 
large), one can always find an elliptical space which fits into the Hubble radius even 
if 1 < f2o + Qao < 1.12. On the other hand, the holonomies of such spaces are Clifford 
translations, so that we can hope to apply the robust PSH method to detect a topological 
signal. This will be the purpose of our subsequent paper. 
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